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ROLE  OF  THE  BOUNDARY  CONDITIONS  IN  THE  PROBLEM  OF 


THE  LINEAR  STABILITY  OF  THE  SEDOV  POINT  BLAST  SOLUTION 
1.  Introduction 

A  highly  localized  release  of  energy  in  an  unbounded  pressureless  ideal  gas 
medium  gives  rise  to  an  expanding  spherical  blast  wave.  Because  of  the  absence 
of  a  characteristic  length  scale  in  the  problem,  the  resulting  flow  is  self- 
similar.  It  consists  of  a  shock  followed  by  a  rarefaction  wave.  If  the  density 
0  of  the  background  medium  satisfies  a  power-law  distribution  (p  ~  R~*  ,  where  R 
is  the  distance  from  the  blast  center),  then  the  solution  of  the  problem  can  be 
reduced  to  quadratures.  This  analysis  ’was  carried  out  by  Sedov  (19^6,  1959)  and 
independently  by  Taylor  (1950)  and  Von  Neumann  (19^+7).  The  analogous  implosion 
problem  was  treated  by  C-uderley  (19^*2). 

Numerous  workers  have  investigated  the  stability  of  the  Sedov  point  blast 
solutions.  In  its  full  generality  the  problem  can  be  stated  as  follows:  liven 
a  particular  3edcv  solution,  specified  by  the  values  of  1  and  the  adiabatic 
index  v,  for  what  wave  numbers  i ,  n  dees  the  deviation  in  the  shock  amplitude 
associated  with  an  infinitesimal  perturbation  about  this  solution  grow  faster 
than  the  shock  radius  5  as  S  -*■  =?  If  no  such  wave  numbers  exist,  the  'unper¬ 
turbed  point  blast  solution  is  stable.  At  present,  a  solution  of  this  general 
problem  is  lacking,  .is  we  shall  note,  very  few  results  of  any  sort  have  been 
established  beyond  iispute. 

Manuscript  approved  Marcn  10.  1986. 


It  is  not  trivial  to  demonstrate  the  stability  even  of  planar  shocks,  a 
problem  treated  originally  by  D'yakov  (195*0  and  later  by  Erpenbeck  (1962).  They 
found  that  for  a  wide  class  of  differing  fluids  (characterized  by  different 
equations  of  state),  planar  shock  waves  propagating  in  a  uniform  medium  are 
stable.  The  physical  mechanism  is  easy  to  describe.  A  small  ripple  in  the  shock 
front  gives  rise  to  divergence  (respectively,  convergence)  in  the  curved  portions 
ahead  of  (behind)  the  main  front.  These  regions  become  weaker  (stronger)  than 
the  unperturbed  shock,  hence  propagate  slower  (faster)  than  average,  thus 
reducing  the  amplitude  of  the  ripple.  Evidently  the  longer  the  perturbation 
wavelength,  the  weaker  the  stabilizing  effect  of  this  mechanism,  so  by  analogy  we 
would  expect  1=0  perturbations  around  the  Sedov  solutions  to  be  the  most 
unstable . 

When  variation  of  the  density  of  the  background  medium  is  taken  into 
account,  however,  this  mechanism  becomes  obscure.  For  2  >  3,  the  unperturbed 
Sedov  shocks  accelerate  with  increasing  radius.  It  is  conceivable  that  the 
ripples  ahead  of  the  shock  can  run  away,  while  those  behind  fall  further  behind, 
leading  to  instability.  Lerche  and  Vasyliunas  (1976)  and  Isen’oerg  (1977),  who 
treated  the  case  y  *  1,  have  claimed  that  for  some  values  of  y,  Sedov  shocks  are 
unstable  for  all  values  of  7.  Their  conclusions  are  remarkable  in  predicting 
instability  at  short  ’wavelength  (i  >>  1).  llevman  (1979)  has  disputed  the 
results,  on  the  ground  that  lerche  and  Vasyliunas  (1976)  and  Isen1  erg  (l?T7) 


improperly  treated  the  boundary  conditions  at  the  she 


cs  i rout • 


Bernstein  and  Book  (i960)  studied  the  stability  of  a  restricted  class  of 


Sedcv  solutions,  those  in  which  the  velocity  behind  the  shock  is  proportional  to 
tne  distance  R  from  the  origin  (homologous  or  uniform  expansion).  This  type  of 


blast  wave  was  apparently  first  considered  by  Primakoff  (see  Ccurant  and 


Friedrions  19^3).  He  found  that  for  0  =  const.,  the  point  blast  similarity 


solution  can  be  written  explicitly  in  closed  form  when  y  =  7. 
generalized  the  results  zo  r.cnuniform  densities  and  found  for 
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the  value  of  3  for  which  explicit  solutions  with  homologous  velocity 
distributions  exist. 

Bernstein  and  3ook  (1980)  linearized  the  fluid  equations  about  the 
Primakoff  blast  solution  and  applied  the  first-order  Rankine-Kugcr.iot  equations 
at  the  shock  location.  The  resulting  eigenvalue  problem  reduces  to  a  quartic 
equation  for  the  stability  index  F,  all  of  whose  roots  have  negative  real  part. 
They  therefore  concluded  that  this  class  of  Sedov  solutions  is  stable, 
contradicting  the  conclusion  of  Isenberg  (1977). 

Subsequently  Gaffet  (I98U)  presented  a  refutation  of  the  work  of  3ernstein 
and  Book  (1980).  Kis  criticism  of  the  latter  focused  on  the  neglect  of  the 
entropy  perturbations  arising  at  the  shock  front.  This,  he  claimed,  resulted 
from  an  incorrect  statement  of  the  perturbed  boundary  conditions. 

To  clarify  the  matter,  we  rederive  the  results  of  3ernstein  and  3ook  (1980) 
taking  into  account  both  the  corrected  boundary  conditions  and  the  entropic 
perturbations.  We  find  that  the  solution  of  the  eigenvalue  problem  is  the  same 
as  that  found  previously,  and  the  conclusion  that  the  modes  are  stable  remains 
unchanged. 

The  plan  of  the  paper  is  as  follows.  We  begin  by  summarizing  the 
followed  by  D'yakov  (195M  in  treating  the  planar  shock  problem, 
because  his  characterization  of  the  acoustic  and  entropic  modes  is  instructive 
the  present  case.  We  then  go  through  the  revised  treatment  of  the  Primakoff 
blast  wave  stability  problem,  using  essentially  the  notation  of  Bernstein  and 
Bock  (198G).  The  paper  concludes  with  a  summary  of  the  results  and  briefly 
explores  their  implications. 
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D'yakov  (195**)  analyzed  the  stability  of  a  planar  shock,  propagating  in  the 
positive  x  direction.  He  chose  as  his  basic  equations  the  adiabatic  law 
expressed  in  terms  of  the  entropy  s,  Euler's  equation,  and  an  equation  for  the 
pressure  p.  The  latter,  taken  together  with  the  adiabatic  law,  is  equivalent  to 
the  continuity  equation.  The  location  of  the  shock  front  is  perturbed  by  an 
amount 

C  =  C  exp[i(kx  -  cut)].  (1) 

o 

All  other  perturbed  quantities  then  are  taken  to  vary  as  exp[i(kx  + 

Zy)i.  In  the  frame  of  the  unperturbed  shock,  the  linearized  equations 
assume  the  form 


(Zv  -  u»)  5  s  =  0  ; 


( Zv  -  u )  >5  v  +  kV<$o  *  0; 
x 


(Zv  -  u)  ov  +  Z'/5o  =  0; 
7 


( Zv  -  <u)  5p  +  oc2(kov  +  Z6 v_  )  =  0, 

x  y 


where  7  =  o "*  is  the  specific  volume,  c  =  (yp/o)l/2  is  the  speed  of  sound,  and 
first-order  fluid  quantities  are  distinguished  with  a  5.  Ercm  Eq.  (2)  it  follows 
that  two  types  of  solutions  are  possible  .indicated  respectively  by  (l)  and  (2);. 

For  type  (1,  Zv  -  u  =  C,  while  5sU  *  0.  It  fellows  from  Eqs.  (3)  and  (<*)  that 
dp*'  =  C,  and  from  Eq.  (5)  that 


/  ,  \  /  ^  \ 
,  ^  [  _  ’  .  '  _  ) 
kc  v  +  z .  ov 

x  X  y 
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"  *,  \r» '  .^rvrv'r  v~»  vnrvw  v>  a  ■%  ^  >  • 


(the  condition  for  incompressibility),  and  the  perturbed  specific  volume 

<SV^^  =  (3V/3s)  5s  is  nonzero.  This  is  the  "entroric"  mode,  which  advects 
P 

entropy  and  vorticity  perturbations  away  from  the  perturbed  shock. 


The  second  type  of  solution  of  Eq.  (2)  satisfies  -  w  *  0,  while  5s 

0.  Multiplying  Eq.  (3)  by  k,  Eq.  (i*)  by  l 2>  adding  and  using  Eq.  (5)  to 
(2)  (2) 

eliminate  5v  '  and  Sv  ,  we  obtain  the  acoustic  disoersion  relation: 
x  / 


(a)  -  i2v)*  =  (k2  +  £2)c2. 


(2)  (2)  (2) 

The  quantities  5v  ,  5v  ,  and  5p  are  then  related  bv  any  two  of  ocs.  i.2)- 
x  i 

(5),  Tor  example 


U9v  -  ui)  5v'2)  +  kV5p'2'  =  C, 


(£2v  -  «)  5v^2)  +  £2V5p(2)  =  C. 


3y  virtue  cf  the  vanishing  of  the  perturbed  entropy  5s  ',  we  have 


5p(2)  *  -(o2/v2)5v(2). 


The  general  solution  of  the  perturbed  equation  consists  of  a  superposition  of 
the  entrocic  and  acoustic  nodes. 


3y  expanding  the  Rar.kine-Kugonict  relations  through  terms  in  first  order, 
four  boundary  conditions  can  be  obtained.  These  are  the  condition  that  the 
tangential  component  be  continuous  across  the  shock, 


two  conditions  obtained  from  the  Jump  in  the  momentum  normal  to  the  shock  and 


■ 


* 


f  ■' 


k 


i.- 


from  the  continuity  equation,  which  D'yakov  (195M  rewrites  in  the  form 


6v^  +  6v 


(2)  _  v  -  v  ('Sp(2)  {V(l)  + 


7 


7 


(12) 


2  p-  p 


v  -  7 


and 


2ifl?  Sp^  <$V(l)  +  <$V(2^ 


(13) 


v  p—  p  V  -V 

and  an  equation  derived  from  the  energy  jump  condition.  Since  for  given  v,  7, 
and  p,  p  is  a  known  function  of  V,  D'yakov  (1956)  takes  this  last  condition  in 
the  form 


5p(2)  =  (3?/37)k  (5V(1)  +  $V(2'), 


(1M 


where  the  subscript  H  stands  fcr  "Kugoniot."  *e  have  consistently  used  an 
overbar  to  indicate  quantities  ahead  of  the  shock. 


Equations  (6)  and  (3)-(lU)  constitute  a  set  of  eight  algebraic  relations  in 


Cl)  (o'*  ( 1  )  fl)  fpl 

the  eight  coefficients  3Vvi,;,  $Vv  ,  $vv  ,  5v;i;,  6v.v<s/. 

x  y  x 


Sv 


(2) 


7 


(o) 

,  5p '  ,  and  s . 


Equation  (T)  is  to  be  understood  as  a  definition  of  Z„  in  terms  of  k  and  2. 


yakov  (195^)  gees  on  to  analyze  the  dispersion  relation 


2ojv  /.  2  .  -o2  >  ,<P  ,  w,  ,  v2  r3V«  . 

-  +  — )  =  ( —  +  )  vo)  -  Z „v )  1 1  +  —  l — ) 

•P  w  3  p 


(15) 


v  ‘r  w 

which  follows,  emphasizing  the  role  of  the  Kugoniot  curve  in  determining  the 


stability  of  the  perturbations.  The  interested  reader  is  referred  to  that  paper 
for  the  details. 


For  the  sake  of  brevity  we  will  not  go  through  the  whole  derivation  and 


solution  of  fluid  equations  linearized  around  the  Primako f f-Sedov  solutions.  In 
this  section  we  will  refer  to  Bernstein  and  Book  (198c)  and  to  equations 
contained  in  it  by  using  the  abbreviation  "3B."  Except  as  noted,  the  notation 
here  is  that  same  as  in  BB. 

The  correct  forms  for  the  Lagrangian  perturbed  density  0.  and  pressure  p  , 

-L  1 

expressed  in  terms  of  the  infinitesinal  displacement  C(R,t)  of  an  element  of 
fluid  whose  unperturbed  trajectory  is  R  (r,t),  are 


0-,  =  P  *  < ) 

and 

?-  =  p(-y  +  X  ), 


:i£) 


(IT) 


where  <  and  X  are  time-independent  quantities  which  were  taken  equal  to  zero  in 
33.  Substituting  (16)  and  (IT)  in  the  perturbed  equation  of  motion  then  yields 


*  r7 


2-v 

r 


V(rV7 


£) 


=  f[< 


(18) 


where  v  =  2  for  the  cylindrical  case  and  v  =  3  for  the  spherical  case.  The 
solution  of  Ec.  (lo)  may  be  written  as  the  sum  of  a  particular  and  a  homogeneous 


where  the  harmonic  dependence  of  the  solution  is  contained  in  H(  )  =  exp(im  i) 


(u=2)  or  H(9,$)  =  (v=3).  Note  that  the  definition  of  3  here  differs 

slightly  froc  that  in  3B. 

Since 

(f3)  =  Sf3-1  ?  +  3  (3-1)  F3'2  f2. 


ve  have  by  virtue  of  33(2.23)  and  (2.2*0 


-2  ,v(y-1)  +  2 


(fP) 


=  -az 


wnere 


3  [2(3-1)  -  v ( Y -1 )  I 
vVY-1) 


writing 


(  )  1 

Xv~;  =  ar  H, 

»• 


r  7 


.  x(s!  =  br3 


we  get  Ec.  33(3.1*0 


(y  *  a)  a  +  (a  +  V  -  1)  -  l] 


( o 
\  •£ 


( o 


(2 


(2 


(2 


free  the  r-ccnponent  of  (it,',  and  an  squatter.  identical  with  3E(3.15/  by  talc  in 


of  (13).  Multiplying  the  first  of  these  'ey  a  *  v  -  1  and 
subtracting  the  second,  ve  obtain  the  equivalent  but  simpler  relation 


A  (—  b  *  a)  +  ,y  -  1;  ((a  +  j  -  1J  a  -  bi  =  0, 


The  "particular"  solutions  of  Eq.  (18),  which  evidently  correspond 
D'yakov's  entropic  mode,  depend  on  time  as  f(t).  Assuming 


ve  can  write 


v  5  “lu 
<  =  K  r  H, 


+  1  r6"1!!. 


?r(P)  ■  *  r5fE, 


r  ?  .  =  3r°fS, 


where  K,  L,  A  and  B  are  constants.  Then  substitution  in  the  r-componer.t  and 
divergence  of  (13)  yields  two  relations  satisfied  by  the  particular  solution. 


(5  +  1)  A  -  3  +  K  =  -  (5  +  v  -  1] 


(1  -  y3) 


*  v  -  1)  [(5  +  1 )A  -  3  +  K]  »  ~  [(6  +  v  -  1)  -  A]  (1  -  oB)  +  AA.  (32) 


Multiplying  (21)  by  <5  +  v  -  1  and  subs-! 


racting,  we  are  ^eft  with 


=  vA  +  y3 


-f  A  *  C;  substitution  in  (32)  then  yields 


K  =  (v  -  2)A  +  3. 


’.■V  A  ^  ^-V  7 


Taking  into  account  the  additional  terms  arising  from  <  and  X  in  Sqs.  (16 
and  (IT),  we  can  rewrite  the  first-order  Hankine-Hugoniot  relations  [3B(3.3^*  - 
(3.36)]  as 


-rcr  +  (v-2)  {r,  -  Zj  + 


ICSJ 


=  qc; 


1  2(;  -  *r) 

:5r  +  v  (y-1)’+  2  Y  +  1  ’ 


2tV. 
ir 


-yra  +  v  (z  -  Z^  +  Xsf 


4  p  S  S  V 


ir 


(y  +  l)  P 


(35 


(36 


(37 


y  -  1 
2 

where  a  =  7  •  ; ,  with  all  expressions  evaluated  at  r  =  s  f  ,  i.e.,  on  the 

shock  front.  If  we  write 


a 


2  s 


(38 


(39 


and  substitute  ;  = 


-  _  ,(?) 


in  las. 


(35)  -  (3' 


we  find 


y  -  1 


^q.  -  >  + 


o  V 


;!sa 


+  3 


[(’J  -  2} a.  +  b 


{va  +  yb  +  (q  -  v )z  - 


2[v(y-1)+21 
Y  +  1 

1 


•}  sQ  fC 


(L  -  vA  -  y3)  s5  f5 


But  Eqs.  (33)  and  (3*0  imply  that  the  right  hand  sides  of  (1»0)  and  (U2)  vanish. 
Thus  all  the  coefficients  associated  with  the  particular  solution  (i.e.,  the 
er.tropic  mode)  cancel  out.  We  are  left  with  three  equations  which,  are  identical 
with  those  found  in  BE.  These  yield  the  boundary  condition 


28a  +  (y  -  1)  b  =  0. 


Equations  (25),  (26),  and  (13)  constitute  three  simultaneous  expressions  for  - 
a/b.  Equating  these,  we  can  solve  for  a,  3*  and  u  as  before,  obtaining  the  same 
quart ic  and  hence  the  same  values  for  the  stability  index 


r  =  [(y  -  1)  «  -  Y  -  1)  +  6.  (WO 

The  arguments  presented  here  must  be  modified  slightly  to  treat  the  case  A  =  0, 
but  the  result  is  the  same  as  that  found  by  taking  the  limit  A  *  0  in  SB (1.8)  and 
(1.13)  which  are  analytic  in  A. 


4.  Discussion 


We  have  seen  that  there  are  indeed  terms  in  the  corrected  analysis  which 
correspond  to  D'yakov's  (195^)  entropic  mode,  hut  they  drop  out  of  the  problem 
and  do  not  alter  our  earlier  result.  Gaffet  (198U)  was  correct  in  pointing  cut 
the  term  corresponding  to  X  *  0  in  Eq.  (17)  for  the  pressure,  but  he  failed  to 
note  the  analogous  term  in  Eq.  (16)  for  the  density.  The  parallel  with  D'yakov's 
(195M  analysis  of  the  planar-shock  (section  2)  is  not  complete;  D'yakov  found  a 
nonvanishing  5V^  (whence  <Sp  ^  *  0),  but  5 =  0.  We  note  that  in  our 
analysis  the  entropic  mode  has  the  same  time  dependence  ~f(t)  as  the  unperturbed 
shock.  Hence  the  density  contrast  <  satisfies  dic/dt  =  0,  whereas  in  the  planar 
case  [by  virtue  of  the  incompressibility  condition,  Eq.  (6)]  it  is  the  perturbed 
density  of  the  entropic  mode  which  satisfies  4—  =  i(Q  -  2-v)  p^  =  0. 

CL  L  c. 

3ernstein  and  Book  (1980)  made  no  use  of  the  tangential  velocity  condition 
analogous  to  Eq.  (11).  In  our  notation  this  can  be  shown  to  take  the  form 

n  x  (  +  7  x  7  ;  =  0,  (4.5) 

where  n  is  the  normal  to  the  unperturbed  shock.  Gaffet  (198M  also  noted  that 
the  perturbation  I  found  by  Bernstein  and  Book  (1980)  dees  not  vanish  at  the 
shock  front,  as  it  should.  We  readily  see  that  Eq.  (45),  together  with 

1  =  0  (1*6) 

at  the  shock  front.  Just  suffice  to  establish  the  connection  between 

a,  b  and  A,  B,  i.e.,  to  determine  the  coefficients  of  the  entropic  part  of  the 

solution  and  the  exponent  5,  which  turns  out  to  satisfy 

5  =  1  +  2s/(y-1 ) .  (1? ) 


The  analysis  presented  by  Bernstein  and  3ook  (1980)  and  here  is  restricted  tc 
the  discrete  modes,  another  point  which  Gaffet  (198U)  criticized.  Although  it  is 
conceivable  that  the  continuum  modes  (which  are  associated  with  the  initial 
conditions)  could  give  rise  to  instability,  we  believe  this  is  unlikely. 

Certainly  most  conventional  problems  of  the  stability  of  ideal  fluids  in  finite 
geometry  can  be  treated  in  terms  of  the  normal  modes  only. 
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